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Abstract 

The computed log(T/y)-log(m/m 0 ) relation for the stellar initial 
mass range, 0.6 < m/niQ < 120 . 0 , and the stellar initial metallicity 
range, 0.0004 < Z < 0.0500, tabulated in an earlier attempt (Portinari 
et al. 1998) is fitted to a good extent by a four-parameter curve, 
expressed by a double power-law, for assigned stellar initial metallicity, 
which can be reduced to a three-parameter curve, expressed by a single 
power-law, for the whole set of stellar initial metallicities. The relative 
errors, i?[log(r/y)] = 1 —log(rfit/y)/log(r/y), do not exceed about 2 % 
and 4%, respectively. The extent to which the interpolation curve, 
expressed by a single power-law, can be extrapolated towards both 
high-mass and low-mass stars, is also investigated. High-mass star 
lifetimes are understimated by a factor less than 2 up to m/mQ = 

1000 and by a fiducial factor less than 4 up to m/m© —)• -|-oo. Low- 
mass star lifetimes are overstimated by a factor of about 3 down to 
m/m© = 0.25 and by an unacceptably large factor down to m/m© = 

0.08. The interpolation curve, expressed by a single power-law, is 
related (in differential form) to a generalization of the equation of the 
exponential decay, which could be the starting point for a theoretical 
interpretation. As a simple application, the star mass fraction of a 
single star generation with stellar initial mass function defined by a 
power-law, is plotted vs. the logarithmic stellar lifetime. The star 
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mass fraction declines in time at a decreasing rate for mild stellar 
initial mass function and at an increasing rate for steep stellar initial 
mass function, where a linear trend is exhibited for a value of the 
exponent close to the Salpeter’s value (pun ~ —2.35). 
keywords - stars: evolution - stars: formation. 
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1 Introduction 


For isolated stars, the stellar lifetime is defined as the time needed to move 
from the zero-age main sequence up to the giant branch and through any 
subsequent giant evolution. There is no easy way to infer stellar lifetime as 
a function of stellar initial mass, from a selected star sample belonging to 
an assigned population. Accordingly, r-m or log(r/y)-log(m/m 0 ) relations, 
where r is the stellar lifetime and m the stellar initial mass, are interpolation 
curves to the results from stellar evolution models (e.g., Maeder and Meynet 
1989; Padovani and Matteucci 1993). For further details, an interested reader 
is addressed to recent overviews (e.g., Romano et al. 2005). 

A systematic dependence of the log(r/y)-log(m/mQ) relation on the stel¬ 
lar initial metallicity, Z, is available in tabular form within the mass range, 
0.6 < m/niQ < 120.0, and the initial metallicity range, 0.0004 < Z < 0.0500, 
where the changes due to diffenent initial metallicities can be neglected to 
a hrst extent (Portinari et al. 1998). Accordingly, lifetime is a strongly de¬ 
creasing function of the initial mass and only weakly dependent on the initial 
metallicity (e.g., Wiersma et al. 2009). 

Interpolating the computed log(r/y)-log(m/m 0 ) relation via a simple, 
continuous, derivable curve, could be useful in problems where the first 
derivatives, dr/ dm or dm/ dr, and higher order derivatives are needed, 
such as chemical evolution models. The current note focuses on a computed 
log(r/y)-log(m/m 0 ) relation where the dependence of stellar lifetime on stel¬ 
lar initial mass and metallicity is exploited (Portinari et al. 1998), aiming to 
the expression of four-parameter curves involving a double power-law, more 
specihcally an exponential whose argument is, in turn, a power. 

To this respect, using standard regression techniques seems to be inap¬ 
propriate, in that different stellar evolution models yield different results for 
an assigned choice of input parameters, due to still persisting uncertainties 
on some physical processes, and standard regression curve procedures would 
be of little signihcance in the case under discussion. Then a different strategy 
must be exploited, where four “crossing” points are specified as belonging to 
the interpolation curve. 

With regard to a hxed stellar initial metallicity, the coordinates of the 
crossing points can coincide with computer outputs, [log(mu/m 0 ), log(ru/y)], 
U = 0, A, B, C, to be hxed as appropriate. With regard to the whole set 
of initial metallicities for which the results are available, the ordinates of 
the crossing points can be determined as average values of their counter¬ 
parts related to the whole set of initial metallicities, while abscissae remain 
unchanged, [log(mu/m 0 ), log(ru/y)], U = 0, A, B, C, to be hxed as appro¬ 
priate, where log(ru/y) = [log(ru/y)i -h ... -h log(ru/y)v]/A^ and A^ = 5 in 
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the case under discussion (Portinari et al. 1998). 

Four-parameter interpolation curves, for both selected initial metallicities 
and the whole set of initial metallicities, are derived in section |2j The results 
are presented in section [3l An illustrative application is shown in section 01 
The discussion and the conclusion are performed in section [5l 

2 Four-parameter interpolation curves 

A complete set of computed stellar lifetime within the initial mass range, 
0.6 < m/iriQ < 120.0, and the initial metallicity range, 0.0004 < Z < 0.0500, 
is available in tabular form (Portinari et al. 1988). More specihcally, stellar 
evolution has been computed for 30 stellar initial masses, where 15 lie within 
the range, 0.6 < m/rriQ < 2.0, 7 within the range, 2.5 < m/rriQ < 9.0, 
6 within the range, 12.0 < m/rriQ < 60.0, with the addition of m/rriQ = 
100.0,120.0; initial metallicities have been hxed to Zi = 0.0004, Z 2 = 0.0040, 
Z 3 = 0.0080, Zi = 0.0200, Z 5 = 0.0500. It is worth noticing Z 4 , initially 
conceived as solar value, is supersolar according to a recent investigation 
yielding Zq = 0.134 (Asplund et al. 2009). 

The computed log(r/y)-log(m/mQ) relation, listed in tabular form in an 
earlier attempt (Portinari et al. 1998), is plotted in Fig. 0] where different 
symbols relate to different initial metallicities: Zi - crosses; Z 2 - diamonds; 
Z 3 - triangles; Z^ - squares; Z^ - saltires. 

The interpolation of the results shall be performed via four-parameter 
curves dehned as (Caimmi 2014): 


y = f{x) = Cl exp{-C 2 x'^) + 03 ] x> 0 ; 


( 1 ) 


where Ci, C 2 , C 3 , 7 , are parameters to be determined by fitting to a good 
extent. 

Regression curves based on standard methods, such as least squares pro¬ 
cedure, could be of little signihcance in that different stellar evolution codes 
with same input parameters yield different results, due to still persisting 
uncertainties on specihc details of stellar evolution theory. 

For this reason, a different strategy shall be adopted, namely interpo¬ 
lation curves are constrained to pass through four crossing points, (xo,|/o), 
{xA,yA), {xByys), {xc,yc), which implies the determination of the unknown 
parameters. More specifically, the substitution of the coordinates in Eq. ([T]) 
yields a system of four transcendental equations in four unknowns, Ci, C 2 , 
C3, 7 , as: 


Vo — Cl + C 3 ; Xq — 0 ; 

yu = Ciexp{-C 2 xl) + C 3 ; U = A, B,C ; 


( 2 ) 

(3) 
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where Xq = 0 , for convenience. 

The combination of Eqs. ([I]) and ([2]) yields: 


y-c^ 

Cl 


Cl - Vo+ y 

Cl 


[exp(-C' 2 )]^^ > 0 ; 


which is equivalent to: 


6*2 


X In 


c-i 


Cl - Vo+ y 


(4) 

(5) 


where Ci < 0 or Ci > |/o — ?/ > 0 in the case under discussion. 

The particularization of Eq. ([5]) to (x{/, yu), (xy, yv), after performing the 
ratio on both sides of related equations, produces: 

pvy \n[{Ci - yo + yw)/Cl] _ 

\xu) \n[{Ci-yo + yij)/Ci\ 

with the above mentioned restrictions for yu, yy, to ensure positive argument 
of related logarithm. In the case under discussion, both x and 7 are real, 
which implies the ratio on the right-hand side of Eq. (jH]) is non negative i.e. 
the arguments of the logarithms are both either larger or lower than unity. 
The related condition reads: 


max(j/u, 2/v) < 2/0 ; min(|/u, 2/v) > 2/o ; (7) 


which holds in the case under discussion, keeping in mind yu < yy < yo with 
no loss of generality, as shown in Fig.[TJ 

Taking the logarithm on both sides of Eq. (| 6 ]) yields: 


. xy .. Cl-yo + yy .. Cl-yo + yu 

7 In — = In In-—-In In-—- 

Xu Cl Cl 


( 8 ) 


where U = A; V = B, C; with no loss of generality. The cross product of 
related equations reads: 


In^ 


, Cl - 2/0 + 2 /b , , Cl - 1/0 + 2 /A 

In In---In In — 


Xk L 
Xk L 


, Xb 

= In — 


Cl 

1 1 Cl — 1/0 + 2 /c ,1 

In In-—-In In 

Cl 


Cl 

Cl — 2/0 + 2 /A 


Cl 


(9) 


which, in particular, can be solved for the crossing points, (0,?/o), (xa,2/a), 
( 1 , 1 / b ), ( 2 , 1 / c ), where xb = 1 , xc = 2 , for convenience. 
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The intersection between the curve on the left and right-hand side of 
Eq. ([H]) yields the value of Ci. The value of can be readily inferred from 
Eq. ([2]) as: 


C3 — Vo — Cl 


( 10 ) 


and the value of C 2 is determined via particularization of Eq. (|3]) to (1,?/b). 
The result is: 


C2 


In 


Cl 

Cl —yo + iiB 


( 11 ) 


hnally, 7 is determined via particularization of Eq. fl 8 |) to U = A: V = B. C. 
The result is: 


7 = 



In In 


Cl-yo + yc 
Cl 


— In In 


Cl —yo + yB 

Cl 


( 12 ) 


where xb/xq = 1/2 in the case under discussion. 

With regard to interpolation curves related to selected stellar initial metal- 
licities, the expression of the mathematical variables, {x,y), in terms of the 
physical variables, [log(m/m 0 ), log(r/y)], is constrained by the domain of the 
interpolation curve, dehned by Eq. (II]), a; > 0 , which implies the following: 


y = log ; (13) 

where mo is the sample star exhibiting the lowest mass i.e. the longest life¬ 
time, mo/mQ = 0.6 in the case under discussion (Portinari et al. 1998). 

The crossing points can be explicitly expressed as: Pq = [0, log(ro. 6 /y)], 

Pa = [-log0.6,log(ri.o/y)], Pb = [1, log(r 6 .o/y)], Pc = [2, log(r 6 o.o/y)], 
where Tm/mQ is the lifetime of a star with selected initial metallicity and 
initial mass equal to m in solar units. 

With regard to interpolation curves related to ordinates of crossing points 
averaged over the whole set of stellar initial metallicity, the expression of the 
mathematical variables, (x, y), in terms of the physical variables, [log(m/mQ), log(r/y)], 
for the abscissa, x, is provided by Eq. (IT^ . while the ordinate, y, reads: 

where Zi, 1 < i < 5, are the stellar initial metallicities for which stellar 
evolution has been computed (Portinari et al. 1998). The crossing points can 


X = 


( m\ I mo 

— I - log ^ 


mp 
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be expli citly expres sed as: Pp = [0, log(ro . 6 /y)], Pa = [-logO. 6 , log(ri,o/y)], 
Pb = [l,log(r 6 .o/y)], Pc = [2,log(r6o.o/y)]. 

Within the stellar initial mass range under consideration, 0.6 < m/mQ < 
120 . 0 , the stellar lifetime range reads 6 < log(r/y) < 11 , which implies ^ 0 
via Eqs. (IT^ and flT^ . Accordingly, use can be made of the relative errors: 

R{y) = l-y^ ; (15) 

y 

where y is inferred from computer output via Eq. flT^ or flT^ . and yfn is the 
counterpart of y on the related interpolation curve. The results are htted to 
a good extent provided relative errors maintain small enough, within a few 
percent say. 

3 Results 

With regard to selected stellar initial metallicities, the interpolated log(m/m 0 )- 
log(r/y) relation is plotted in Fig.[T] as full curves, which exhibit negligi¬ 
ble difference with the exception of the supersolar stellar initial metallicity, 

Z = Z 5 , where slightly lower values are shown. A dotted horizontal line, 
log(r/y) = 8 , marks a conventional boundary between long-lived (up) and 
short-lived (down) stars, which could be of interest for simple chemical evo¬ 
lution models (e.g., Pagel and Patchett 1975; Caimmi 2011) and type la 
supernova progenitor models (e.g., Acharova et ah 2013). 

The parameters of the interpolated log(m/mQ)-log(r/y) relation, Ci, C 2 , 
( 7 * 3 , 7 , are listed in Tabled! The parameters lie within narrow ranges, namely 
4.92 < Cl < 5.24; 0.78 < C 2 < 1.00; 5.49 < C 3 < 5.71; 0.96 < 7 < 1.06. 
The special case, 7 = 1, makes the interpolation curve reduce from a double 
power-law (an exponential whose argument is, in turn, a power) to a simple 
exponential. The interpolation curve related to averaged ordinates of crossing 
points and its reduction to a simple exponential, cases 6 and 7 of Table dl 
are plotted in Fig. [1] as a dashed and a dotted curve, respectively. 

The relative errors, R{y), dehned by Eq. fd5l) . as a function of the loga¬ 
rithmic stellar initial mass, log(m/mQ), are plotted in Fig.[2|for interpolation 
curves related to both selected stellar initial metallicities, cases 1-5 of Table 
dl (top panel) and averaged ordinates of crossing points, cases 6-7 of Table 
dl (bottom panel), where different symbols are captioned as in Fig.dl An 
inspection of Fig.[2| shows the relative errors do not exceed about 2% in the 
former alternative and less than 4% in the latter one. Then it may safely be 
thought interpolation curves, expressed by Eq. ([T]), with values of parame¬ 
ters listed in Tabled! satisfactorily £t to the computed log(r/y)-log(m/m 0 ) 
relation. 
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Table 1: Parameters (p) of the interpolated log(m/m 0 )-log(r/y) relation, Ci, 
1 ^ 2 ; C's; 7; fo^ Selected stellar initial metallicities, Z = Z 1 -Z 5 with regard to 
cases (c) 1-5, and for ordinates of crossing points averaged over the whole set 
of stellar initial metallicity, case 6 , with the addition of the special choice, 
7 = 1, leaving the remaining parameters unchanged, case 7. See text for 
further details. 


P: 

c 

Z 

Cl 

C 2 

^^3 

7 

1 


5.10551 

0.785151 

5.52594 

0.96950 

2 

^2 

5.23694 

0.784218 

5.49141 

0.96686 

3 

^3 

5.20510 

0.824914 

5.60580 

0.98915 

4 

^4 

5.08428 

0.904558 

5.81444 

1.03125 

5 


4.92130 

0.995420 

5.93483 

1.05302 

6 

z 

5.07899 

0.862584 

5.70612 

1.00322 

7 

z 

5.07899 

0.862584 

5.70612 

1 


The special case, 7 = 1, deserves further attention in that the interpola¬ 
tion curve reduces to a simple exponential via Eq. ([T]), as: 


?/= Ciexp(-C2a;) + C3 ; 

and the substitution of Eq. flT^ into flTBD after some algebra yields: 


r 


m 




log r = — 


+ C?, 


J, 


VflQ 


C[ = Cl exp 


C2log(^ 


^0 

m© 


= 0.6 ; 


2 In 10 


(16) 

(17) 

(18) 


which implies stellar lifetime is expressed by a double power-law, more specif¬ 
ically an exponential whose argument is, in turn, a power, as: 


r 


y 


= expio 




m 




-c' 


+ C 3 


(19) 


where the parameter values are: 

C( = 4.19439 ; C 2 = 0.374616 ; C 3 = 5.70612 ; (20) 


according to case 7 listed in Table [H 










4 A simple application 


With regard to a single star generation, let the stellar initial mass function 
be a power-law as: 


0 




( 21 ) 


where A is the star formation efficiency and —3 < p < —2 for the cases of 
interest, in particular p = —2.35 relates to a classical investigation on the 
solar neighbourhood (Salpeter 1955). 

By dehnition, diV = (j){m/mQ) d{m/mQ) is the number of stars born 
within the mass range, m/rriQ^ d(m/m 0 )/ 2 , and the initial mass of the star 
generation reads: 


M*(0) = AniQ 


myif/rriQ ^ 
nn^tlmQ 7710 



( 22 ) 


where 7T7Mf and ?T 7 jnf are the upper and lower stellar initial mass limit, re¬ 
spectively, in solar units. 

The mass of the star generation (without stellar remnants) at the time, 
t = T^m/mQ,), reads: 


M* 



Amp 


0 


pmImQ 


m 

7770 



(23) 


where T{m/mQ) is the lifetime of a star of mass m in solar units. 

The fractional mass of the star generation, s(r/y) = M*(r/y)/M*(0), 
can be determined from Eqs. fl 22 l) and fl25D after performing integration. The 
result is: 


^ /r\ _ (777/7710)^+^ - ( 7 ? 7 mf/ 7 T 70 )^+P 

Vy/ (n7Mf/n70)2+P - (777n,f/7T70)2+P 

/r\ ln(m/m^f) _ „ = _2 ■ 

\y) ln(777Mf/n7mf) 


p^-2 ; (24a) 

(24b) 


regardless of the star formation efficiency. 

The fractional mass of a star generation, expressed by Eq. fl2TD . is plotted 
in Fig. [3] for p = —3.0, —2.9,..., —2.0, from top to bottom, with regard to 
the interpolation curve expressed by a simple exponential, Eq. (ITB]) . case 7 
of Table [H Related stellar lifetimes, corresponding to the whole range for 
which stellar evolution has been computed in the parent paper (Portinari 
et ah 1998), are marked by dotted vertical lines where mmi/mQ = 0.6 and 
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myii/rriQ = 120.0. A restricted number of stellar initial masses are captioned 
on the top of the box. 

An inspection of Fig. [3] shows the expected trend: mild stellar initial mass 
functions {p ~ —2) imply a larger fraction of short-lived stars and, in turn, 
a change in star mass fraction at a decreasing rate. On the other hand, 
steep stellar initial mass functions {p ~ —3) imply a larger fraction of long- 
lived stars and, in turn, a change in star mass fraction at an increasing rate. 
Interestingly, the transition case where the s(r/y)-log(r/y) relation is linear, 
is close to the Salpeter’s value, p = pun ~ —2.35. 

To gain more insight, the above mentioned special case relates to a 
straight line passing through the points, [log(rmf/y), 0] and [log(rMf/y), 1], 
which can be expressed as: 


T 


log(r/y) - log(rmf/y) 


jJ log(rMf/y) - log(rmf/y) 


(25) 


where = r(minf), XMf = T(mMf), for brevity. The combination of Eqs. fl24ap 
and fl25|) . after little algebra, yields: 



2+p 


log(rMf/y) - log(rmf/y) / m \ 

(mMf/m0)2+P - (mmf/mQ)2+P ym© j 
_ log(T'Mf/y) - log(rmf/y) / 


2+p 



; (26) 


where p = pun via Eq. fl25ll . 

The comparison of Eq. 
of the following relations: 


(mMf/m0)2+P - (mmf/m0)2+P y m© / 

with fliyp term by term implies the validity 




-2-p ; 

log(rMf/y) - log(rmf/y) 


C = 

^ (mMf/m0)2+p - {m^i/mQY+P 


C, = 


log(rMf/y) - log(rmf/y) 


{niMi/mQy+P - {m^i/mQY+P \ mQ ) 


+ log 


\ y 


(27) 

(28) 

(29) 


where, in the case under discussion, the left-hand side of Eqs. fl7rp - fl29p is 
known via Eq. fl20p and rUmf = 0.6, rriMi = 120.0. Then p = pun can be 
evaluated via Eqs. fl20l) . ([27]), and Eqs. fl28p and fl2^ are merely to be verified. 
The result is: 

piin = -2.374616 ; (30) 

where the values of C[ and coincide with their counterparts expressed by 
Eq. fl20p . as expected. 
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In conclusion, a star generation where the initial mass function is ex¬ 
pressed by a power-law within the range, 0.6 < m/rriQ < 120 . 0 , implies a 
linear dependence of the star mass fraction, s(r/y), on the logarithmic stel¬ 
lar lifetime, log(r/y), for an exponent of the initial mass function close to 
Salpeter’s value, piin ~ —2.35. 

5 Discussion and conclusion 

Computed stellar log(r/y)-log(m/m 0 ) relations are interpolated to a good 
extent by double power-laws and, in particular, by a single power-law, Eqs. ([T]) 
and ([TB]) . respectively, regardless of the initial stellar metallicity, within the 
stellar initial mass range, 0.6 < m/niQ < 120.0. The case of a single power- 
law is expecially attractive due to its simplicity over an extended mass do¬ 
main. 

By comparison, a log(r/y)-log(m/m 0 ) relation can be used for massive 
{m/rriQ > 6 . 6 ) stars (e.g., Matteucci and Greggio 1986), while a cumbersome 
expression is necessary for less massive stars (e.g., Renzini and Buzzoni 1986). 
In alternative, four different log(r/y)-log(m/m 0 ) relations can be used within 
the mass range, 1.3 < m/mQ < 120.0 (Maeder and Meynet 1989). A three- 
parameter (depending on the stellar initial metallicity) £t to earlier results 
from the Padua group (Along! et ah 1993; Bressan et ah 1993; Bertelli et 
ah 1994) is restricted to Z < 0.03 and the extrapolation does not reproduce 
the correct trend as m/mQ —>• -|-oo (Raiteri et al. 1996). 

With regard to the interpolation curve, expressed by Eq. flT 6 |l . the follow¬ 
ing relations hold for stellar initial masses outside the domain of computed 
stellar evolution (Portinari et al. 1998) used in the current attempt: 




( 31 ) 


m 


( 32 ) 


- y -f-cx) ; 

rriQ 


which is the correct trend, in the sense that (failed) stars with mass, m/mQ < 
m-iai/mQ ~ 0.08, never start hydrogen burning and related lifetime can be 
conceived as inhnite. On the other hand, stars with increasingly high mass ex¬ 
hibit hnite lifetime regardless of initial metallicity. The comparison between 
computed lifetimes (Portinari et al. 1998) and their counterparts inferred 
from the interpolation curve, expressed by Eqs. flTBl) . fl2U]) . case 7 of Tabled], 
is shown in Table [2] It is apparent lifetimes calculated via the single power- 
law interpolation curve outside the domain towards high stellar masses are 
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Table 2: Comparison between computed stellar lifetimes (Portinari et al. 
1998), Tpa 98 , and their counterparts inferred from the single power-law inter¬ 
polation curve, case 7 of Table [H tj, outside the domain towards high-mass 
stars. See text for further details. 


m/rriQ 

'TPags/y 

T-y/y 

150 

3.0E6 

2.2286E6 

200 

2.7E6 

1.9162E6 

300 

2.5E6 

1.5894E6 

500 

2.0E6 

1.3032E6 

1000 

2.0E6 

1.0507E6 

-|-oo 


5.7061E5 


understimated with respect to computed values, by a factor not exceeding 2 
up to m/rriQ = 1000 . 

Concerning stars with decreasingly low mass, a similar comparison is 
shown in Table [31 where computed stellar lifetimes are from a different source 
(Laughlin et al. 1997), restricted to stellar initial metallicity, Z = 0.02 
(Laughlin and Bodenheimer 1993). It is apparent lifetimes calculated via 
the single power-law interpolation curve outside the domain towards low 
stellar masses are unacceptably overstimated with respect to computed values 
but, interestingly, the discrepancy is within the same order of magnitude for 
m/rriQ = 0.25, where stellar evolution still attains the giant phase (Laughlin 
et al. 1997). Accordingly, an absence of the giant phase would imply a 
different log(r/y)-log(m/m 0 ) relation. 

A simple but acceptable linear £t to the results listed in Table |3l leaving 
aside m/niQ = 0.25, reads: 

which implies log(ri.o/y) = 10.800 and log(ro. 6 /y) = 11.257, larger but of 
about the same order of magnitude with respect to their counterparts deter¬ 
mined from stellar evolution (Portinari et al. 1998). 

The computed log(r/y)-log(m/m 0 ) relation, represented in Fig.Hl ex¬ 
tended to the results listed in Tables [ 2 | and [31 is shown in Fig. [31 where, 
in addition, the single power-law interpolation curve and the interpolation 
line, Eqs. flTB]) . f[2U]) . and ([33]) . respectively, are also plotted together with 
the asymptotic limit of the interpolation curve for inhnite stellar initial 
mass. The global trend exhibited by the computed log(r/y)-log(m/m 0 ) re- 
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Table 3: Comparison between computed stellar lifetimes (Laughlin et al. 
1997), rLa 97 , and their counterparts inferred from the single power-law inter¬ 
polation curve, case 7 of Table [H ry, outside the domain towards low-mass 
stars. See text for further details. 


m/iriQ 

'TLagr/y 

A/y 

0.08 

1.1E13 

3.24E16 

0.10 

6.2E12 

4.40E15 

0.12 

4.2E12 

9.72E14 

0.14 

3.4E12 

2.93E14 

0.16 

2.7E12 

1.10E14 

0.18 

2.0E12 

4.78E13 

0.20 

1.8E12 

2.35E13 

0.25 

1.7E12 

5.71E12 


lation is reminiscent of earlier results restricted to stellar initial metallicity, 
Z = 0.02 (Maeder and Meynet 1989), where the transition between two dif¬ 
ferent regimes occurs above (e.g., Wiersma et al. 2009) instead of below 
(Fig.Sj) the solar mass. 

In conclusion, the single power-law interpolation curve, expressed by 
Eqs.®, ([20]), can safely be extrapolated towards high-mass stars (m/m© > 
120 . 0 ) yielding understimated values within a factor of two up to m/m© = 
1000 and within a hducial factor of four up to m/m© —>• -|-cx3. Addi¬ 
tional caution must be used for the extrapolation towards low-mass stars 
(m/m© < 0 . 6 ), yielding overstimated values within a factor of about three 
down to m/m© = 0.25 and to unacceptably larger overstimates down to 
m/m© = 0.08. 

The question if the single power-law interpolation curve, expressed by 
Eqs. fllOp . 0201 ) . arises from a mere coincidence or has a theoretical interpre¬ 
tation, is outside the aim of the current attempt but some considerations 
can be performed. Differentiating both sides of Eq. fllOp . after little algebra 
yields: 

^ = -A-A-'*/ ; (34) 

T 777 

i = - ■ rfi = — A = C'(C'' ; /3 = C'+1 ; (35) 

y mQ 

where the special case, (3 = 0, reduces to the equation of the exponential 
decay, / being related to the number of decaying nuclides and rfi to the time. 
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Accordingly, Eq. fl5^ learns (i) longer stellar lifetimes imply larger change of 
stellar lifetime with stellar initial mass and vice versa, and (ii) higher stellar 
initial masses imply lower change of stellar lifetime with stellar initial mass 
and vice versa. In the author’s opinion, a theoretical explanation of Eq. ffTHjl 
should necessarily start from the above mentioned points. 

The fractional stellar mass as a function of the logarithmic stellar lifetime, 
shown in Fig. [21 relates to a single star generation and to a power-law stellar 
initial mass function, which are limiting cases. In reality, multiple star gen¬ 
erations appear even in low-mass stellar systems and the stellar initial mass 
function exhibits a more complex trend with respect to a simple power-law. 
On the other hand, the ideal situation depicted in Fig.|3]is expected to show 
a similar trend with respect to the real case, at least for environments where 
the whole amount of pristine gas was turned into stars in a short time, less 
than 1 Gyr say, such as the inner halo, the thick disk, and globular clusters. 
It is worth noticing the stellar mass is less than the total mass, in that the 
contribution of stellar remnants (white dwarfs, neutron stars, black holes) 
and baryonic dark matter is neglected: in any case, the related contribution 
is expected to be of a few percent at most. 

With regard to globular clusters with age equal to ro.g (hrst vertical line 
on the right of log(r/y) = 10), an inspection of Fig. [3] shows a present star 
mass fraction, 0.07 ~ s ~ 0.35, according if —2 > p > —3, in particular 
s Ri 0.15 for the Salpeter’s exponent, p = —2.35, close to a linear dependence 
of the star mass fraction on the logarithmic stellar lifetime. Then a mild 
stellar initial mass function, p ^ — 2 . 1 , implies an initial mass larger by a 
factor of about ten, for globular clusters with age equal to ro.g, in the case 
under consideration. 

The main results of the current note can be summarized as follows. 

(1) The computed log(r/y)-log(m/m 0 ) relation for the stellar initial mass 
range, 0.6 < m/niQ < 120 . 0 , and the stellar initial metallicity range, 
0.0004 < Z < 0.0500, tabulated in an earlier attempt (Portinari et 
al. 1998) is htted to a good extent by a double power-law for assigned 
initial metallicity, which can be reduced to a single power-law for the 
whole set of initial metallicities. The relative errors, i?[log(r/y)] = 
1 — log(rfit/y)/ log(r/y), do not exceed about 2 % for the double power- 
law and 4% for the single power-law. 

(2) The interpolation curve, expressed by a single power-law, can be ex¬ 
trapolated towards both high-mass and low-mass stars. In the former 
alternative, stellar lifetimes are understimated by a factor less than 2 
up to m/rriQ = 1000 (Portinari et al. 1998) and by a hducial factor less 
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than 4 up to m/niQ —)■ +oo. In the latter alternative, stellar lifetimes 
are overstimated by a factor less than about 3 down to m/mQ = 0.25 
and by an unacceptably large factor down to m/mQ = 0.08 (Laughlin 
et al. 1997). 

(3) The hrst derivative, d(r/y)/d(m/m 0 ), related to the single power-law 
interpolation curve, is a generalization of the equation of the exponen¬ 
tial decay, which could be the starting point for a theoretical interpre¬ 
tation of the log(r/y)-log(m/m 0 ) relation. 

(4) In the special case of a single star generation with stellar initial mass 
function dehned by a power-law, using the single power-law interpola¬ 
tion curve within the mass range, 0.6 < m/mQ < 120 . 0 , the star mass 
fraction declines in time at a decreasing rate for mild stellar initial mass 
function {p ~ — 2 ) and at an increasing rate for steep stellar initial mass 
function [p ~ —3), where a linear trend is exhibited for a value of the 
exponent close to the Salpeter’s value (pun ~ —2.35). 
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Figure 1: Computed log(r/y)-log(m/m 0 ) relation for 30 different stellar ini¬ 
tial masses within the range, 0.6 < m/mQ < 120.0, in connection with 5 
different stellar initial metallicities (Portinari et al. 1998). Different sym¬ 
bols relate to different metallicities as indicated. Four-parameter interpola¬ 
tion curves are shown as full for selected stellar initial metallicity; dashed 
for ordinates of crossing points averaged over the whole set of stellar initial 
metallicity; dotted for parameters, Ci, C 2 , C 3 , as in the last case and 7 set 
equal to unity. The dotted horizontal line marks a conventional boundary 
between long-lived (top) and short-lived (bottom) stars. See text for further 
details. 
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Figure 2: The relative error, Riy) = 1 — log(rfit/y)/ log(r/y), as a function of 
the logarithmic stellar initial mass, log(m/mQ), with regard to interpolation 
curves related to selected stellar initial metallicity (top panel) and to averaged 
ordinates of crossing points (bottom panel). Different symbols correspond 
to different stellar initial metallicity as shown in Fig.[T] See text for further 
details. 
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Figure 3: Fractional mass of a star generation, s(r/y) = M*(r/y)/M*(0), 
as a function of the logarithmic stellar lifetime, log(r/y), for a power-law 
stellar initial mass function with exponent equal to (from top to bottom) 
p = —3.0, —2.9,—2.0, with regard to a computed log(r/y)-log(m/mQ) 
relation interpolated by a simple exponential, restricted to the domain, 0.6 < 
m/mQ < 120.0. Related stellar lifetimes, corresponding to the whole range 
for which stellar evolution has been computed (Portinari et ah 1998), are 
marked by dotted vertical lines, a restricted number of which are captioned 
on the top of the box. 
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Figure 4: Computed log(r/y)-log(m/m 0 ) relation, plotted in Fig.[Tl extended 
to both lower (Laughlin et al. 1997) and larger (Portinari et al. 1998) stellar 
initial masses. Lower masses relate to stellar initial metallicities, Z = 0.02. 
An interpolation straight line and power-law curve, respectively, are also 
shown (dotted). The asymptotic limit of the interpolation curve towards 
inhnite stellar initial mass is marked by a dotted horizontal line. See text for 
further details. 
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